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^ , Abstract 

G . Phase-field simulations of the martensitic transformation (MT) in the austenitic 

'^ \ matrix, which has already undergone the plastic deformation, are carried out. For 

W ■ this purpose the elasto-plastic phase-field approach of incoherent MT developed in 

the previous work [Kundin et. al. J. Mech. and Phys. Solids 59 (2011) 2012] is used. 

The evolution equation for the dislocation density field is extended by taking into 

account the thermal and athermal annihilation of the dislocations in the austenitic 
K* . matrix and the athermal annihilation at the transformation front. It is shown that 

the plastic deformation in the austenite caused by the MT interacts with the pre- 



cn 

Q"^ ' deformed plastic strain that leads to the inhomogeneous increasing of the total 

t^^ ■ dislocation density. During the phase transformation one part of the dislocations in 

^" , the pre-deformed austenite is inherited by the martensitic phase and this inheritance 

depends on the crystallography of MT. An other part of dislocations annihilates 
at the transformation front and decreases the dislocation density in the growing 
martensite. Based on the simulation results a phenomenological dependency of the 
inherited dislocations on the martensitic fraction and the plastic deformation in the 
martensite and austenitic matrix is proposed. 
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1 Introduction 



The morphology of the martensitic microstructures is strongly related to the 
mechanical properties of steel and alloys. It is important to predict the for- 
mation of martensitic microstructure precisely. Furthermore, during the MT 
the incremental elastic strain energy is reduced by the formation of a hetero- 
geneous array of different orientation variants of the martensitic phase and by 
plastic deformation (plastic accommodation). The development of appropriate 
models allows to calculate the strains associated with the phase transforma- 
tion. The understanding and control of localized plasticity and its interaction 
with the transformation front is quite essential e.g. for TRIP-steels and shape 
memory alloys [T|2|3] . 

Recently, phase-field (PF) models have been extensively studied as a powerful 
tool for predicting microstructural evolution and applied to the martensitic 
transformation. Khachaturyan and co-workers developed the phase-field mi- 
croelasticity (PFM) theory pi5|6] . which integrates microelasticity into the 
phase-field model employing the fast Fourier transform algorithm. The model 
has further applied to investigate the MT in single crystals [7] and poly- 
crystalline systems [8] as well as for the MT with external stresses. It has 
been successfully applied to various coherent phase transformations includ- 
ing the prediction of many complicated strain-induced morphological patterns 



Besides the coherent phase transformations the MT in technical alloys is usu- 
ally associated with plastic strains. The effect of plasticity, which can be rep- 
resented as the generation and the motion of dislocations on the MT, is very 
complex. On the one hand, if the MT is caused by applied loads, the evolu- 
tion of the plastic deformation reduces the driving force of the transformation. 
On the other hand, a local plastic relaxation allows the accommodation of the 
strain caused by the MT and the nucleation and growth of the new martensitic 
variants. Finally, statistically stored dislocations are an irreversible structural 
change which affect the energy landscape by there own elastic fields. Thus, a 
general phase-field model for MT should include not only the driving forces 
originating from the elastic fields, but also the effects resulting from the plastic 
strain fields. 

The plastic activity in the phase-field theory has already been treated by mod- 
eling the dislocations [H] and their dynamics P^16|17|18|19] . A phase-field 
model of the evolution of a dislocation system based on the time depending 
Ginzburg-Landau (TDGL) equation was developed by Wang et al. [17]. At 
the same time, Koslowski et al. [19] formulated a phase-field theory of the 
dislocation dynamics for an arbitrary number and arrangement of dislocation 
lines based on energy minimization. 



The plastic deformation can also be added to the phase-field model by in- 
troducing a plastic strain field defined at the mesoscale. A version of this 
approach has been recently proposed by Zhou et al. pUII^ where the plastic 
strain is related to the inter-dislocation distance. The evolution of the disloca- 
tion phase fields is described by the TDGL equation similar to the evolution 
of the martensitic phase fields. The difference to the individual dislocation 
dynamic models is that the model proposed by Zhou et al. describes the evo- 
lution of a dislocation mixture. Another version of this approach, which take 
into account the strain energy relaxation due to plastic accommodation, has 
been proposed by Yamanaka [2^P5] based on the PFM theory and the elasto- 
plastic PF model suggested by Guo et al. [2^ . 



To describe the plastic deformation in single crystals many models use the 
crystal plasticity framework in which the evolution of plastic strain is described 
by means of an elastic driving force. It is the so called phenomenological visco- 
plastic approach [25|26|27|28|29j . Crystal plasticity methods developed by Ma 
and Roters [30|25|29j have recently made remarkable progress and allow the 
numerical study of deformation processes on the basis of thermally activated 
dislocation evolution. The model preserves the crystallographic features of dis- 
location slip processes and captures the commonly accepted concepts of dis- 
location processes in plastic deformation, especially various dislocation inter- 
action processes as interactions of mobile and statistically stored dislocation, 
the formation of locks and dipoles or the thermal and athermal annihilation in 
a continuum dislocation density framework. The original concept [30] for fee 
single crystals has been extended to polycrystals, considering grain boundary 
interactions and geometrically necessary dislocations as well as the extension 
to dislocations in bcc crystals [3T|32|33|29|25] . It should be noted that a num- 
ber of other works have been published recently which couple the phase-field 
simulation with crystal plasticity schemes at the mesoscale [Mf55] and do not 
take into account the plastic accommodation caused directly by the solid-solid 
transformation. 

Crystal plasticity models, their further developments based on microscopically 
interpretable state-variables and their evolution based on sound physical prin- 
ciples are very attractive for coupling with solid-solid phase transitions and 
offer new research directions and new insights into the mechanisms of marten- 
sitic phase transformations. One contribution in this direction is the work of 
Kundin et al. [36], where the crystal plasticity model [30|25|29] is coupled 
to the phase-field model of the MT that allows to resolve the dynamics of 
individual martensitic plates. It has been demonstrated that the formation 
of special martensite morphologies (butterfiy type) in Fe-Ni is assisted by the 
formation of plastic strain fields in the austenitic matrix, which then are inher- 
ited by the growing martensitic plates as it was experimentally observed in the 
works [37]. The model also uses methods of earlier models for the simulation 
of transformation induced plasticity 



The main goal of the current paper is the understanding of the evolution of the 
dislocation during the MT and the interaction of the plastic strain with the 
transformation front. The kinetic of the MT is described by the elasto-plastic 
PF model of the MT. The dislocations are treated as a separate dislocation 
field and their evolution is treated using the model of thermally activated 
dislocation motion [30|25|29] . In our investigation we concentrate on the in- 
heritance of the dislocations by the martensitic phase. The previous work in 
this direction is the paper of Ostwald et al. |10], where the dislocation inheri- 
tance is considered as a result of the interaction between phase-transformation 
and plasticity effects. The authors introduced a so-called plastic inheritance 
law, which describes the plastic strain evolution during the transformation 
front propagation. Is was assumed that the dislocations neither generate nor 
annihilate on the transformation front, i.e. the overall amount of plastic defor- 
mations remains constant. The dislocations can be only inherited in a growing 
martensitic phase or pushed by the transformation front, increasing the plastic 
strain in the parent austenitic phase. An inheritance probability function was 
suggested for the quantitative description of the inheritance process, which is 
equal to the fraction of the inherited dislocation density. In contrast to this 
inheritance mechanism we proceed from the fact that there are no physical 
driving force to push dislocations from the transformation front. We con- 
sider the mechanism of the annihilation of the dislocations on the propagating 
martensitic front and the generation of the dislocations in the austenitic ma- 
trix. The probability of the dislocation inheritance is defined as the fraction 
of the remaining dislocation density after the annihilation. 

We first start with recalling the phase-field microelasticity model of the MT 
extended by taking into account plastic strains in Section 2. The differential 
equation describing the evolution of plastic strain is presented in Sections 2.3. 
The coupling of phase-transformation and plasticity effects is incorporated in 
one algorithm. Details on the numerical investigation of the transformation 
front and the plastic strain are provided in Section 3, where the model is 
applied to a hypothetical Fe-Ni alloy. In addition, when various phase-field 
model parameters are applied, the resulting effect on the inheritance kinetics 
was predicted by the suggested model and compared with the reported data. 



2 Theoretical part 

2. 1 The main parameters of the elasto-plastic phase-field model of MT 



Following the concept developed in the PFM theory [5f7f6llT3] we consider a 
coherent multi-phase mixture with the order parameters r]p{r) and the local 
stress- free strain tensor £°j(p, r) = e'^j{p)rjp{r) where p identifies a crystallo- 



graphic variant. The order parameters are treated as phase- field variables so 
that rjp G [0, 1] and satisfy the property I]p=i^p(i") = 1 ~ 0a(i') where 0a(i") 
is the phase-field variable responsible for the presence of the austenite. 

The dislocations are characterized by a Burgers vector b*^") and a slip plane 
with a normal n^"^ defined by the crystallography of the austenitic (or marten- 
sitic) phase. Here, a is an index corresponding to a slip system. We use the 
symmetric Schmid tensor for a slip system a defined as M*^") = n^"-' <S) ni^°'^ 
where m^"") = h^°') /b^""^ is the unit vector in the direction of the Burgers vector, 
which expresses the slip direction. Then we define a dimensionless dislocation 
function through the dislocation densities of mobile, pM,a, and immobile, pi^a, 
dislocations 



r) = 6WJp,,„(r)+pM,„(r) (1) 



as an analog of a phase-field variable. It is equal to the unity if the distance 
between the dislocations in the full volume is equal to the Burgers vector. The 
local stress-free plastic strain tensor caused by dislocations of an a-th slip 
system is given then by 



4(a,r) = M,,(a)0^(r). (2) 

Note that the sum of all dislocation functions is smaller than 1 X!a=i </'a(r) ^ 1. 
We assume that z/ is the number of the martensitic variants and /i is the 
number of the slip systems. 



2.2 The evolution of the order parameters in the elasto-plastic phase-field 
model 



The formulation of the total free energy of a system as a function of the order 
parameters is the key step in the development of phase-field models. In the 
standard phase-field formulation the total free energy is divided into three 
terms: the gradient energy term, the double well potential term (both terms 
are responsible for the interface energy) and the chemical free energy that pro- 
vides the chemical driving force for the transformation. Concerning the MT, 
there are an additional elastic strain energy that suppresses or accelerates the 
transformation. The chemical free energy and the double well potential term 
in PFM theory of the MT are usually considered together and approximated 
by a Landau polynomial expansion with respect to the order parameters 



In the present study we reformulate the PFM theory in terms of the standard 
phase-field modeling [52|53] . The evolution of the order parameters in the 
system with the plastic strain is governed by the elasto-plastic phase-field 
kinetic equation 



1_^ = KV\{r) - Hf'^^ + A/,; + ^-^^'4.(^)4 



CijkiSiAp) 



f: elM (r/,(r) - r/,) + f: Mh(/3) {^{v) - 0^ 



Lg=l /3=1 



Ciifc«£L(p)e*^m,(e)<TL(k)ene'''-"— -^, (3) 



(2vr) 
where the Fourier transform of the stress is calculated by 



ikr 73 



d'r. 



<^L(k) = / c^n.k E ^]M Mv) -%) + Y. MjkiP) (0"(r) - 4 

•'^ U=l /3=1 

(4) 
Here Cijki is the elastic modulus, Qjm is the Fourier transform of the Green 

function, f/ and (p are the volume averaging functions, afj'^ is the applied 

stress. 

The first term in the right hand site of eq. ([3]) is the gradient term, which 
forces interfaces to have a finite width. The second term is the double potential 
function and the third term is the driving force. Then A^ is the mobility of the 
interface, ii" = ^ is the gradient coefficient, if = -^ is the coefficient before 
the double well function, ^ and 7 are the interface width and the interface 
energy, respectively, which are assumed to be equal for all solid-solid interfaces. 

The functions /' and g' are the derivatives of the model functions with 
respect to the chosen order parameter. The model function / is the double 
well potential and is defined as follows 



fiv} = Ev',-2j:vt+{j:vU ■ (5) 

9=1 9=1 V9=l / 

The function g is responsible for the chemical driving force and defined as 




^(^1=4E^'-3 E^ll • (6) 

9=1 \Q=1 

These model functions serve the goal to keep the order parameters between 
and 1. In Appendix A we show that the chosen model functions satisfy the 
properties of the standard phase-field model and derive their relation to the 
Landau polynomial expansion usually used in PFM theory. 

The undercooling responsible for the chemical driving force is given by 

Af = QM{TM-T)/TM, (7) 
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where Qm and Tm are is the latent heat and T is the temperature of the MT, 
respectively. 



2. 3 The evolution of the plastic strain 



In this section we describe the time evolution of the dislocation density gener- 
ated during the MT and forming plastic deformation strain. The dislocation 
density will first evolve in the austenitic phase. Then, due to the interaction 
of the plastic strain with the transformation front, one part of the dislocations 
will cross the interface without or with small changes of the slip system and 
be inherited in the martensite. The other part of the dislocation will annihi- 
late by the interaction with the dislocations of the same slip system collected 
at the transformation front. The kinetics of these processes depends on the 
growth velocity of the transformation front and the resistance stress of the 
matrix. It should also be taken into account that a part of the dislocations 
can annihilate in the bulk phase. In the following we consider a single crystal 
system, i.e. a system which contains only one grain without boundaries. 

In order to simulate the evolution of the dislocation density and thereby the 
dislocation function, we define additionally to the mobile and immobile dislo- 
cations the parallel dislocation density p-p^a and the forest dislocation density 
Pf,« for a slip system a following [SD] as 



2A;bT 

PM,a — T;;7Ta/PF,oPP,o, 

CiC2C2,LrO^ 



pp,a = E Pi,/3 1 Sin (n("), n(/^) x m^^)) | , (9) 

/3=1 



= 1 



PF,a = E Pi,P |cos (n("), n(« X m(«) | , (10) 



where G is the shear modulus, /cb is the Boltzmann constant, T is the tem- 
perature and Ci,C2,C3 are material constants. 

Following the theory presented in [25|30] we calculate the evolution of the 
immobile dislocation density as 



Pl,a = C4y/pF^-fa - C5pl,ala " C7 exp i^-^j^ ) YT^^^'""' ^° 

- C9P/,„ E ^P + Clo0AV[na X V]p/,a, (11) 



where 7^ is the plastic shear rate, Qbuik is the activation energy for chmb and 
Tr, is the external stress. 



la 



Eq. flTTl) combines three processes in the bulk phase: immobilization of the 
mobile dislocations, non-thermal annihilation with the immobile dislocations 
of the same slip system and thermal annihilation by climb of edge dislocations 
with material constants C4, C5 and cy, respectively. The probability of these 
processes are governed by the plastic shear rate 7q,. 

Furthermore, we propose to take into account an additional process of the non- 
thermal annihilation at the martensitic transformation front with a constant 
cg. The probability of this process is proportional to the local dislocation den- 
sity and the velocity of the MT front, which is related to the evolution of the 
order parameter 77^. Physical meaning of this term is that the transformation 
front collect the dislocations, which can later annihilate with the dislocations 
of the same slip system during the front propagation. The increasing collection 
or striking of the dislocations on the front should increase with the increasing 
transformation rate, because the system has less time to relax. In the following 
we will call the annihilation at the MT front as the annihilation term. In this 
study the two cases are considered with and without the annihilation term. 
Note that the parameter cg is a phenomenological parameter, which should 
be found experimentally. In the simulation it is is chosen of five time smaller 
than the parameter C5, which is responsible for the annihilation in the bulk 
phase. 

Similar to the gradient energy formulation used in the phase-field model of 
individual dislocations [17] , the gradient of the dislocation density pi^a is added 
to the elastic energy of the dislocation and subsequently the Laplace term is 
added to the evolution equation (ITT]) . This gradient can be attributed to the 
core energy of the dislocations. The coefficient Cio has to guarantee a smooth 
transition of the deformation strain field profile across the sheared non-sheared 
martensite-austenite boundaries. The value of cio should be the same order as 
the square of the dislocation core size. A similar term was used in work [21] 
for the simulation of the 7'-coarsening in the presence of the dislocation field 
and is reasonable. 

The plastic shear rate can be defined according to the Orowan equation as 

7« = 4>APM,abaVa, (12) 

where we assume that the plastic deformation is hindered in the martensite 
due to much higher fiow stresses compared to the austenite. We describe the 
evolution of plastic slip during the process only in the austenitic phase and 
multiply the plastic shear rate with (pA- 

The average dislocation velocity f „ based on the model of thermally activated 



dislocation motion can be found as [30|25f29] 

Va = A„i/„exp ( --^j exp f -^^j^Vaj , for |r"| > r^^,,, (13) 

where Tp^sg = ciGha^/p-p,a + PM,a is the passing stress of the mobile dislo- 
cations, z/q is the attack frequency, Q^^p = Gh\/2 is the effective activation 
energy for the dislocation slip, Aq, = -p^= is the jump width and Va = cs^^Aq, 
is the activation volume, where Ci, C2 and C3 are material constants. 

In the simulations numerical problems arise with eq. (fT3!l due to the very 
strong exponential function of the stress. For numerical reasons we prefer an 
alternative variant of the evolution equation according to a Norton type flow 
rule 

Va = Kv. exp (^-^ j (^-^ j ' (14) 

where Tcut = 0.2 GPa and n = 5 for cubic systems. This type of equation was 
used in [M] to study rafting in Ni-based superalloys. 

The external stress r", which is the driving force of the dislocation evolution, 
is caused by the applied force, growing martensitic phases and dislocation 
fields. It can be derived from the total elastic energy as 



ra = -CijkiMiAa 






+ / Q,HMM(a)e,fi,„,(e)a^Jk)e„e(*'^-'-)— -^ + atfMi^{a), (15) 

JV ylTiy 

where the Fourier transform of the stress (j^„(k) is defined by (j4]). 

The set of equations ([S]) - (IT^ . dSD, (fT5]) are calculated in each step of the 
simulation for all slip systems, then the dislocation functions are defined ac- 
cording to eq. ([1]) and inserted in the kinetic equation ([3]). 



3 Numerical simulation 



3.1 Choice of parameters 



For the numerical investigation we chose a hypothetical Fe-30%Ni alloy, where 
the occurring lenticular-type martensite is characterized by the mixture of 
twinning and dislocation structure. 



Table 1 

Material and processing parameters used in the simulation 



Parameter 


Value 


Temperature of MT 


Tm = 400 K 


Interface energy 


7 = 1.9 X 10-2 J m-2 [51] 


Latent heat 


Qm = 3.5 X 10^ J m-3 U] 


Pre-deforming 


pj = 1012 m-2 


External applied stress 


alf' = 4^/ = -1 GPa 


Shear modulus 


G = 28 GPa 


Poisson coefficient 


i^ = 0.374 


Material constants 


ci = 0.18 




C2 = 5.0 




C3 = 5.0 




C4 = 8.0 • 10^ m-i 




C5 = 10 




C7 = 1 X 10^ m^s'^s 




C8 = 0.3 




eg = 2 




cio = 0.1/^ 


Attack frequency 


uo = 10^0 s-1 


Activation energy 


Qsiip = 2.3 X 10-19 J 


Activation energy 


Qbuik = 2.4 X 10-19 J 


Lattice constant 


ao = 3.59 X 10-1° m 



The material parameters and phase-field model parameters used in the simu- 
lations are listed in Tables 1, 2. 



The reason of the choice of the discretization size is that we first chose the 
relation between the interface width and the maximal capillary length do = 
7/(ai-Eo) = 1-3 X IQ-ii m. That implies that ^/do < 200, which is a require- 
ment of the asymptotic limit in the phase-field model [S2]- In the simulations 
we used the phase-field equations in dimensionless form by measuring length 
in units /q, time in units of tq ( it can be estimated from the experimental 
formation rate for the MT of order 10-^ s) and the energy in units of Eq, 
which is the typical strain energy of MT. Then the mobility of the interface is 
calculated as A^ = 1/{tEq). 
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Table 2 

Phase-field model parameters used in the simulation 



Parameter 


Value 


Units 


Length scale, Iq 


1.3 X 10^9 


m 


Time scale, tq 


3.3 X 10-1° 


s 


Energy scale, Eq 


3.07 X 10^ 


Jm-3 


Discretization size. Ax 


1 


h 


Time step. At 


0.125 


To 


Interface width, ^ 


1.5 


lo 


Gradient coefficient, K 


0.0152 


Eol'o 


Double-well potential coefficient, H 


0.0067 


Eo 


Landau energy parameter, 02 


0.0134 


Eo 


Interface mobility, Ai 


1 


iroEor' 


Undercooling, A/ 


0.06 


Eo 



For these parameters the Landau energy parameter 02 = 27/(01^) in the 
chemical energy function (see Appendix A) equals to 02 = 0.0134 i?o- The es- 
timation of the gradient coefficient in eq. (^ yields K = 7^/(ai) = 0.0152 EqI^. 
For the comparison with work [13] the gradient coefficient is in the same or- 
der [K = 0.01624 ii^o^o) but the Landau energy parameter is ten times larger 
(02 = 0.312 ii^o)- The interface width used in the simulation was very much 
smaller than the chosen discretization size C, - 



IK/H = 0.32/0. This means 
that the increasing parameter 02 results in a decreasing interface width. In 
the simulation we used the parameters listed in Table 1. For comparison we 
also carried out the simulations with 02 = 0.312 E'q (■C = 0.32 /g) to show the 
influence of this parameter on the annihilation term. 



For three variants of the MT we used the Bain transformation matrices: 

£11 = 0.1322, £22 = 0.1322, £33 = -0.1994, variant 1; 
eu = -0.1994, £22 = 0.1322, €33 = 0.1322, variant 2; 
£11 = 0.1322, £22 = -0.1994, £33 = 0.1322, variant 3. 



In the fee lattice of the austenitic matrix 24 variants of the dislocation slip 
were chosen with slip planes of type {111}^ and slip directions of type (101)a- 
For the calculation of the Burgers vector we used b = ao/\/2, where Oq is the 
lattice constant. 
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3.2 Numerical investigation of the dislocation inheritance 



In this section we present the simulation results for the case where disloca- 
tions exist in the pre- deformed austenite. A nucleus containing 3 martensitic 
variants is generated in a system of size 38 x 38 x 38 Ax and 64 x 64 x 64 Ax, 
respectively. Periodic boundary conditions are applied to the systems for all 
fields. 

An external force of 1 GPa is imposed in direction [101]. This force suppresses 
one martensitic variant and promotes two other variants (variants 1 and 2), 
which produce a twinning structure. The dislocation density evolves during 
the simulation due to the applied stress and the stress from the MT. 

The probability of dislocation inheritance from the austenitic matrix to the 
martensitic phase is defined as a ratio between the averaged dislocation density 
with and without the annihilation term in eq. ( ITTj) 

T> u\ ^^'° + Pi,adt 



Pi, a + Pl,adt + Cgpl^aVdt' 



where rj = J2p=i % ^^^ ^ identify a space averaged variable over the interface 
between martensitic and austenitic phases 



/o^''°" xr](j)AdV 



X = '\ "'^'"^"'^ . (17) 



Thus, we consider only the dislocations which can take part in the annihilation 
process at the MT front. Note that from the probability Va{t) the contribution 
of the annihilation term to the total dislocation density is derived as 1 — 'Pa(t). 

During the simulation we additionally control the ratio between the mean 
dislocation densities in the martensite (all variants) and austenite 



n^{t) = ^. (18) 

Pl,a\A 

The simulated microstructure during the MT at 1500, 2500 and 5000 time 
steps is shown in Figs. 1-3 for three cases. In the first and second tests (Figs. 
1 and 2) we simulated the microstructure with the interface width $, = 1.5 /q 
in boxes of size 38 Ax and 64 Ax. The microstructure consists of twinned 
martensitic plates, which belong to one martensitic lath. The size of the system 
influences the number and the thickness of the growing plates. It can be seen 
that for the smaller system size the number of the martensitic plates is larger 
and their thickness is evenly distributed in the volume. For the smaller system 
size the thickness of martensitic plates becomes 2-3 times larger at the end of 
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the growth. In the third test (Fig. 3) the interface width was chosen ^ = 0.32 Iq. 
Furthermore the coordinate system was rotated to [112], [111], [110]. In the 
simulation box of size 64 Ax 6 nuclei of the martensitic phase were inserted 
randomly with the mean time interval equal to 100 time steps. One nucleus 
grows in one martensitic lath, which consist of the thin twinned martensitic 
plates. The simulation test shows that for the small interface width ^ = 0.32 Iq 
the MT can not be completed at the chosen undercooling. To reach the full 
martensitic transformation the simulation was carried out at the increasing 
undercooling from A/ to 3A/ for 2500 simulation time steps that corresponds 
to the cooling rate 2 ■ 10^ K/s. 

The 2D sections (x — y-directions) of the microstructure and the dislocation 
density field (in logarithmic coordinates) are presented in Figs. 4 and 5 for the 
first and second test, respectively. It can be seen that the dislocation density 
increases first near the transformation front in the earlier stages of MT and 
strongly increases in the full austenitic volume at the end of the MT. The 
annihilation decreases the dislocation density in the martensitic phase (last 
column in Figs. 4,5). 

For the analysis we chose two slip systems as examples of low and high in- 
teraction stresses between the martensite and the dislocation field. As a char- 
acteristic parameter of this interaction a maximum interaction stress can be 
proposed defined as 

Tinter(a,p) = CijkiMij{a)eli{p), (19) 

which influences the kinetics of the MT and the evolution of the dislocation 
field due to the phase transformation (see eq. US])- It is obviously that the inter- 
action stress between the martensitic plates and the dislocation field depends 
on the orientation relationship between the corresponding martensitic vari- 
ants and the dislocation slip systems. The full list of the interaction stresses 
between 12 slip systems and two martensitic variants is presented in Table 3. 
As an example for the visualization and the qualitative analysis we chose two 
typical slip systems which have different summary interaction stresses (see last 
column of Table 3). It is also remarkable that the evolution of the dislocations 
is faster in the martensitic variant 2 as in variant 1, because the interaction 
stress is not zero for both chosen slip systems. In Figs. 4,5 the difference be- 
tween the diffusion fields in the areas of both martensitic variants can be seen, 
which is large in the variant 2 (yellow color). 

The evolution of the mean dislocation density during the simulation is plotted 
in Fig. 6 (a,b) for two slip systems. It can be observed that the dislocation 
density in the martensite increases stronger without than with annihilation. 
In is obviously that the annihilation decreases the dislocation density in the 
martensite according to eq. (ITT]) . For the slip system 1 the dislocation density 
is larger due to the larger summary interaction stress. 
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Table 3 

Interaction stress between the dislocation slip systems and mertensitic variants. The 



ms used in the simulation are unc 


er lined. 




Slip system 


Mart, variant 1 


Mart, variant 2 


sum 


(iii)[no] (1) 






0.0736807 
0.0736807 


0.0736807 
0.0736807 


(111)[110] 


(111)[011] 


-0.0736807 





-0.0736807 


(111)[011] 


-0.0736807 





-0.0736807 


(lll)[10l] 


0.0736807 


-0.0736807 





(lll)[10l] 


0.0736807 


-0.0736807 





(111)[I10] 





-0.0736807 


-0.0736807 


(lll)[110] 





0.0736807 


0.0736807 


(111)[101] (2) 


-0.0736807 
0.0736807 


0.0736807 
-0.0736807 






(111)[101] 


(1I1)[0II] 


0.0736807 





0.0736807 


(111)[011] 


-0.0736807 





-0.0736807 



The time evolution of the probability of the inheritance, Va-, is plotted in 
Fig. 7 (a). A comparison with Fig. 6 shows that the larger the dislocation 
density and their evolution rate, the larger the probability Va- The simulation 
with the small interface width ^ = 0.32 Zq shows that the annihilation of the 
dislocations at the transformation front is very small. Indeed in this case the 
phase-field model gives maximal 2 discretization points on the interface, while 
for ^ = 1.5 /q the phase- field model gives minimal 4 discretization points. It 
is obviously that for a thicker interface the annihilation processes are more 
intensive. 

The probability Va as function of the martensitic fraction is plotted in Fig. 
7 (b). It slightly decreases for a long time and then strongly decreases at the 
middle of the MT due to the dislocation evolution in the austenite caused by 
the MT. At the end of the transformation the probability strongly increases 
to 1 for the system size 38 Ax. That means that the annihilation of the dislo- 
cations at the interface increases with the increasing dislocation density and 
decreases with the decreasing transformation rate. For the system size 64 Ax 
Va is smaller than for the system size 38 Ax and it does not reach 1 at the end 
of the transformation. It can be explained by the fact that for the larger sys- 
tem size and correspondingly larger distance between the martensitic nucleus 
the content of the retained austenite is larger and the dislocation density is 
smaller (see Fig 6). 
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In Fig. 7 it can be also seen that for the shp system 1 the average probabihty 
of the inheritance is higher than for system 2. This can be explained by the 
higher interaction stress that causes the intensive dislocation density evolution 
and the corresponding decrease of the contribution of the annihilation term. 

In the work [IQ] a convex inheritance probability function was proposed, where 
it was required to match the inheritance probability function to at x = 
and to 1 at a; = 1. It was also assumed that if x tends to zero, the dislocation 
density in the austenite takes high values and forces to be inherited by the 
MT front. To proof the proposed convex inheritance probability function we 
found that the curves in Fig. 7 (b) can be approximated in their increasing 
part by a function of the martensitic fraction, x, similar to the proposed 



Vo, = ko + {l-h)xe^Y>(—Y^)^ (20) 

where ko is a minimal probability and ki is a hardening coefficient. For x = 1 
the probability becomes 1. The value kix is responsible for the increasing 
plastic strain in the retained austenite due to the MT. For the simulated 
curves of the system size 38 Ax the parameters of eq. fl20l) are the following: 
ko = 0.9710, 0.9490 and ki = 0.012, 0.010 for the slip systems 1 and 2, 
respectively. For the system size 64 Ax the dependencies in Fig. 7 (b) can not 
be approximated by the chosen function. 

The time dependency of the ratio between the dislocation densities in the 
martensite and austenite is plotted in Fig. 8 for two slip systems. The behavior 
is similar to the probability of the inheritance. With increasing dislocation 
density in the austenite the ratio TZ decreases. The annihilation term decreases 
the dislocation density in the martensite at the earlier stages and does not 
influence the further evolution. It is interesting that without annihilation the 
dislocation density is higher in the martensite and the ratio TZ is higher than 
1 all the time during the MT. 

The comparison of the kinetics of the MT with and without dislocation evo- 
lution and with and without annihilation term is shown in Fig. 9. The gen- 
eration of the dislocation during the MT increases the rate of the MT and 
the final martensitic fraction. The annihilation on the contrary decreases the 
transformation rate. It can be seen that for small system size 38Ax the final 
martensitic fraction also decreases with the annihilation term. 

From the simulation results it follows that the contribution of the annihilation 
at the transformation front is smaller than 6%. That means that the contri- 
bution of the annihilation in the total process of the evolution of the plastic 
strain should by very small. But in spite of this the influence of the annihi- 
lation on the MT is essential. It is expressed in a decreasing rate of the MT 
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and a decreasing final martensitic fraction. And in addition tlie anniliilation 
decreases tlie dislocation density in the martensitic phase almost to half. 



4 Conclusions 



In the present study the interaction of the MT front with the plastic strain 
in a Fe-Ni alloy has been investigated in detail by means of the elasto-plastic 
phase-field model. 

The model takes into account the formation of accommodation dislocations 
in the austenitic matrix and the annihilation of the dislocations at the grow- 
ing martensitic transformation front. It resolves the effects of the dislocation 
dynamics on the kinetics of the MT. It is found that the part of the disloca- 
tions inherited in the martensite can be defined as functions of the martensitic 
fraction and the plastic strain or the dislocation density in the austenitic ma- 
trix. As result of the numerical simulations the following phenomenological 
dependencies for the probability of the inheritance of the dislocations by the 
martensite are found: 

• 1) Va depends on the plastic strain, the plastic shear rate and the transfor- 
mation rate; 

• 2) Va depends on the orientation relationship between the martensitic vari- 
ants and dislocation sleep systems; 

• 3) Va depends on the mean distance between the nuclei of the martensitic 
phase; 

• 4) Va can be describe as a function of the martensitic fraction. 

Moreover, from our simulation follows the dependency of the probability of 
the dislocation evolution and the inheritance on the orientation relationship 
between the martensitic variants and dislocation slip systems. 



5 Acknowledgement 

Appendix A: The chemical energy function 

The model function / is defined as follows 



q=\ q=l \q=l 
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and its derivative is equal to 



V 



4 = 2hp-3r/^ + 2r/pX]r/2|. (.2) 

/ reduces in the case of one martensitic variant to the standard double well 
function 



f = r)\\-rif (.3) 

with 

/; = 2r/(l-r/)(l-2r/). (.4) 

For this function the numeric constant a\ = \/2/3 [53] . 
The model function g is defined as 



2N 

9iv) = 12 \ It Vt-\(±f,] I (.5) 

with a derivative 



gl = 12Ll-VpT.^v'l (-6) 

Further g reduces in the case of the one martensitic variant to the model 
function 



with 



^; = 127/2(1 -7/). 



It can be seen that the second and third terms in the kinetic eq. ([3]) build 
together the well known Landau-type polynomial 
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h = Hf{ff)-Afg{ff) 

= \<^^il< - i(3«2 + 12A/)X:r/J + \ {2a, + 12A/) \±rf^ | (.9) 

^ 9=1 -^ q=l ^ \q=l 

where the Landau energy parameter a, = 2H = ^ is related to the nucleation 
barrier. 
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Fig. .2. 3D images of simulated microstructures dm'ing MT for the box of size 64 Ax 
with ^ = 1.5^0 at 2500 (a) and 5000 (b) time steps; the surface view at 5000 time 
steps (c). 
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Fig. .3. 3D images of simulated microstructures during MT for the box of size 64 Ax 
with S, = 0.32 Iq at 2500 (a) and 5000 (b) time steps; the surface view at 5000 time 
steps (c). 
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Fig. .4. Simulated microstructure and dislocation field evolution during MT in the 
2D cross-sections of a 3D box of size 64 Ax along [001] axis with ^ = 1.5 ^o- First 
column contains the microstructures at 1500, 2500, 3500 and 5000 time steps, red 
and yellow areas represent the 1 and 2 martensitic variants. The corresponding 
dislocation density field in logarithmic scale is shown in the second column without 
annihilation and in the third column with annihilation term. 
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Fig. .5. Simulated microstructure and dislocation field evolution during isothermal 
MT in the 2D cross-sections of a 3D box of size 64 Ax along [001] axis with S, = 0.32 ^o 
through 6 nucleus. First column represents the microstructure at 1500, 2500, 3500 
and 5000 time steps, red and yellow areas represent the 1 and 2 martensitic variants. 
The corresponding dislocation density field is shown in the second column without 
annihilation and in the third column with annihilation term. 
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Fig. .6. Time evolution of the dislocation density function for the slip system 1 (a) 
and 2 (b). 
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Fig. .7. Time evolution of the probability of the dilocation inheritance for slip sys- 
tems 1 (a) and 2 (b). 
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Fig. .8. Ratio between the dislocation densities in the martensite and austenite as 
function of the martensitic fraction with and without annihilation and diffusion for 
slip systems 1 (a) and 2 (b). 
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Fig. .9. Evolution of the martensitic fraction for two variants (a) and (b) with and 
without annihilation. 
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